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Abstract 

o, 

We present a new staggered discretization of the Dirac operator. 
| In comparison with standard staggered fermions, real and imaginary 

parts are located in different nodes. Doubling gives only a doublet 
of Dirac fermions which we propose to interpret as a physical (lepton 
or quark) doublet. Contrary to usual staggered fermions, we have no 
pi 1 exact chiral symmetry but obtain a generalization of the Ginsparg- 

O . Wilson relation. 



1 Introduction 



In this paper we present a new discretization of the Dirac equation. In com- 
parison with staggered fermions M it creates not four but only two flavours 
of Dirac fermions. This has been reached by placing not only different spin 
components, but also their real and imaginary parts into different nodes. 
These sixteen steps of freedom (two fermions) can be understood, in some 
sense, as the result of "doubling" of a real scalar step of freedom ip(n) on the 
lattice. 

Moreover, these two fermions live on different sub-lattices (ijj on "odd" 
nodes, ip e on "even" nodes), thus, we obtain also a single fermion (eight steps 
of freedom) on the lattice by omitting half of the lattice nodes. But, if the 
neutrino is a Dirac particle, all fermions appear in doublets of Dirac particles. 
In this context, the appearance of a fermion doublet in this discretization may 
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be not a bug but a feature which allows to explain why Dirac particles appear 
in such doublets. 

In our approach the complex structure is an operator among others. 
Moreover, there is no natural complex structure, but, instead, a quaternionic 
structure. The nature of this structure is a topic of future research. 

Once a complex structure is defined, we have to define the operator 7 5 
on the grid. Now, there is no exact chiral symmetry. Instead, we obtain 
formulas for some operators 7 5 on the grid which define a generalization of 
the Ginsparg- Wilson relation. 



2 A Real Representation Of The Dirac Alge- 
bra 

We forget - for some time - about the complex structure. Instead of the usual 
representations with four complex fields, we use an eight-dimensional real 
representation of the operators 7 M defined here by their linear combination 
with d u : 



7 °a -7^ = 



do 


0i 


d 2 




d 3 








di 


-d 




d 2 




d 3 






d 2 




-do 


"01 






d 3 






-d 2 


d l 


do 








d 3 


d 3 








-do 


-0i 


-d 2 






-d 3 






0i 


do 




-0 2 






-d 3 




d 2 




do 


0i 








-d 3 




d 2 


-0i 


-0o 



def 



V 

In the context of this representation, it seems also natural to define (by 
their linear combination with scalar parameters rrii) the following operators 
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The following operator equation holds: 



(7°<9o - f dt + m^f = -□ + 5 ij m imj (3) 

This can be easily seen - this operator iterates three times, in each coor- 
dinate direction, the same trick:Q 

( {m - a,), im '-A ]I ) 2 - ^ + <"•• + a -)<"'- - 9 -« ( I °, ) W 

As follows immediately, the 7^ define a representation of the Dirac ma- 
trices, and the matrices f3 l fulfil the following anticommutation relations: 

+ (3 j f3 i = 5 ij (5) 

and anticommute with all 7^: 

/?V + 7"/? = (6) 
It is also easy to see (and to generalize to arbitrary dimension) that 



7 °(7 1 /3 1 )(7 2 /3 2 )(7 3 /3 3 ) = l. (7) 

3 Complex Structures 

In this representation we have not yet defined a complex structure. There 
are several things in the standard approach to Dirac fermions which depend 
on them. 

First, in the standard approach a Hermitean scalar product (., .) is widely 
used. We have only a standard Euclidean scalar product (., .) yet. Now, for 
a complex structure i these notions are closely related in a simple way: 

The Hermitean scalar product defines an Euclidean scalar product by 

(^0) = i((?M> + (0,V>> (8) 

lr This observation also suggests how to iterate this construction to arbitrary dimension. 
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so that (iijj,i(fi) = ("0,0). For a complex structure i, i 2 = — 1, with this 
property this Hermitean scalar product is defined by the Euclidean scalar 
product by defines the Hermitean scalar product 

<V,0> = (V>,0)- (9) 

Thus, we should not care about the Hermitean scalar product, the Eu- 
clidean scalar product is all we need. Note also that for a complex lin- 
ear operator A, which is a real linear operator with [A,i] = 0, the Her- 
mitean adjoint operator A + and the Euclidean adjoint operator A* cooin- 
side: (A*i/j,(f)) = (ip,A(f)). As a consequence, the classical properties of the 
7-matrices 

= yy/ (io) 

are equivalent to 

(7T = 7°/7° (11) 

These properties are fulfilled in our representation for the standard Eu- 
clidean scalar product (., .) in M 8 . 

Next, the "classical" operator 7 s = — Z7°7 1 7 2 7 3 also depends on the com- 
plex structure. A natural replacement which does not depend on it - the 
expression 7°7 1 7 2 7 3 - we denote with 1: 

t =de f 7V7V = P W lY + ^l = (0 2 = -l (12) 

But, of course, we have to introduce a complex structure if we want to 
connect the Dirac fermion in the usual way with gauge fields. The properties 
we need for a complex structure i are i~ x = i* = —i and [7^, i] = 0. Now, an 
interesting point is that there are several candidates for such a structure: 

i =(3 l (3 2 =l(3 3 (13) 

j =P 2 P 3 =lP (14) 
k =/3 3 /? 1 =LfJ 2 (15) 

which together define a quaternionic structure :0 

2 The classical representation ij = k can be obtained using reverse signs for i,j, k, but 
we prefer this sign convention because it gives 7 s = /3 3 . 
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ij = —ji = —A;; jk = —kj = —i; ki = —ik = —j; i 2 = j 2 = k 2 = — 1 (16) 

For each candidate % for a complex structure, we obtain an own operator 
7 5 —def —it- Especially for i = i(3 3 we obtain 7 s = —i/3 3 i = j3 3 . Thus, it 
seems that to fix the complex structure we somehow have to break spatial 
symmetry, prefer one spatial direction. 



4 Discretization Of The Dirac Equation 

Our representation is appropriate for a discretization of the Dirac equation 
on a regular hyper-cubic lattice. It can be obtained in a quite simple way: 
We start with a naive central difference approximation 



dii)(n) -»• — (i/;(n + ai) -if>(n- o<)). (17) 

This naive discretization leads to the problem of "fermion doubling" . The 
continuous limit of this set of discrete equations gives not only the original 
Dirac equation, but also additional, highly oscillating components, the "dou- 
blers" . In classical computations such doublers may be often ignored, but in 
quantum computations, where the number of steps of freedom is important 
(Pauli principle) this is no longer possible. We obtain in each direction a fac- 
tor to, thus, 2 4 = 16 doublers. Fortunately, eight pairs of doublers decouple 
in a really simple way: It is sufficient to hold only one of the eight real com- 
ponents if) a per node. On the three-dimensional reference cube (£1,62, £3), 
Si G {0, 1} we obtain the following locations for the eight components: 



1J) located at (0,0,0), 

if) 1 located at (1,0,0), 

if) 2 located at (0, 1,0), 

ip 3 located at (1, 1, 0), 



if) 4 located at (0,0,1) 
if; 5 located at (1, 0, 1) 
if; 6 located at (0, 1, 1) 
if) 7 located at (1, 1, 1) 



What remains are sixteen steps of freedom (eight steps of freedom on 
two time steps which we need because of our use of central differences in 
time) which corresponds to a doublet of Dirac fermionsf]. Note that our dis- 
cretization may be interpreted as a way to discretize the d'Alembert equation 

3 This is a variant of a well-known approach to solve the "fermion doubling" problem 
- so-called staggered fermions [|r). The standard staggered grid approach reduces the 
doublers only by a factor four. Because we ignore the complex structure of the standard 
representation, we are free to place "real" and "imaginary" part of the complex fields into 
different nodes. This gives the additional reduction by factor two. 
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for a single scalar step of freedom (p(n) with central differences, which gives 
2 4 = 16 doublers. 

Now, the last doublet decouples too, but in a slightly less trivial way: We 
can distinguish "even" and "odd" nodes on the full space-time lattice. The 
central difference equations on even (odd) nodes connects only values on odd 
(even) nodes. Thus, we obtain two fermions ip e and ip a on even resp. odd 
nodes. On the four- dimensional reference cube (eq, ei, £2^3), £% £ {0,1} we 
have 



ifj e ° located at (0,0,0,0), 
if) e l located at (1, 1, 0, 0), 
ifj e 2 located at (1, 0, 1, 0), 
if) e 3 located at (0, 1, 1, 0), 



-0 e 4 located at (1, 0, 0, 1) 
-0 e 5 located at (0, 1, 0, 1) 
V> e 6 located at (0,0,1,1) 
i/j e 7 located at (1, 1, 1, 1) 



ifj o located at (1, 0, 0, 0), ^o 4 located at (0, 0, 0, 1) 

ip l located at (0, 1, 0, 0), if) 5 located at (1, 1, 0, 1) 

i) 2 located at (0, 0, 1, 0), ^ 6 located at (1, 0, 1, 1) 

ip 3 located at (1, 1, 1, 0), i/j 7 located at (0, 1, 1, 1) 

But, instead of removing one sub-mesh to describe a single Dirac fermion, 
we propose to accept above doublers as a way to describe a physically mean- 
ingful flavour doublet. Remarkably, if the neutrino is a standard Dirac par- 
ticle, then all fermions of the standard model appear in doublets. The ap- 
pearance of a fermion doublet in our approach may be, therefore, not a bug 
but a feature which allows to explain the existence of these doublets. 

4.1 Lagrange Formalism 

To define the Lagrange formalism on the lattice note that the equations for 
the ip are located at the nodes of i\) e and reverse. Therefore, we can use the 
Lagrangian 



(20) 



L = £>eWo) = -X>We) (21) 

and obtain the equation for Dip / e = as the Euler-Lagrange equation 
for ipe/o- We can also rewrite this Lagrangian as 

£=^5>(-ir 7 o (^) (22) 
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4.2 Fermion Families and Lattice Distortions 

We have observed that to fix the complex structure we somehow have to 
break spatial symmetry. But there is a simple way out of this. Instead 
of one scalar step of freedom <p(n) on the lattice, we can consider a vector 
field - thus, three components Lp l (n). Now, each component has a natural 
"preferred direction" and, therefore, a natural complex structure. 

Moreover, a vector on a lattice is a quite natural step of freedom. It is, for 
example, the natural way to describe lattice distortions with a shift vector 
field u l (ri). 

On the other hand, in the standard model we have three fermion families 
- three copies of each fermion. This suggests to explain on the kinematic 
level the three fermion families using the hypothesis that the fundamental 
steps of freedom of the "theory of everything" are three-dimensional vector 
fields u l {n) of spatial distortions. 

5 Chiral Symmetry On The Lattice 

One problem with standard staggered fermions is that they have the wrong 
number of doublers (four) to allow a natural physical interpretation in the 
standard model. In our approach, we have only a pair of Dirac fermions, 
and pairs of fermions appear in the standard model as quark pairs as well as 
lepton pairs (if the neutrino is a Dirac particle). 

The other problem is that there is exact chiral 7 5 symmetry on the lattice. 
As a consequence, the doublers appear in pairs with reverse chiral charge. 
This does not fit the situation in the standard model (cf. ||). Now, in 
our approach we do not have exact chiral 7 5 symmetry. Instead, we have 
a replacement for this symmetry. This replacement fulfills properties which 
define a generalization of the famous Ginsparg- Wilson (GW) relation ||. 

Let's assume now that one of the complex structures, namely i = l/3 3 , 
has been chosen. To understand chiral symmetry we have to define the op- 
erator 7 5 = [3 3 on the lattice. Note that it cannot be anymore a pointwise 
operator as for Wilson fermions and staggered fermions - it connects com- 
ponents which are located in different points. Now, we propose to consider 
the following operator as a candidate for 7 s on the lattice: 

(7 5 0)(n ei)en ) = (f){n even - h z ), (7 5 0)(™odd) = <\>{n odd + h z ) (23) 

It is easy to see that it approximates the continuous 7 s . More interesting 
is that some exact properties remain valid: 
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(7 5 )* = 7 5 , (7 5 ) 2 = 1 
We can also define, as an alternative, the operator 7 s by 



(24) 



(7 s '(p)(n even ) = (j){n even + h z ), (j 5 <P)(n odd ) = (p(n odd - h z ) (25) 
Similarly, we obtain 

(f7 = f\ (7 5 ) 2 = 1 (26) 

If we define the operators V, O by 7 5 = rfV = 7 5 (1 — aO) 
we obtain the Ginsparg- Wilson (GW) relation for O: 

7 5 + 7 5 = aO^O (27) 

Moreover, we have also the following important commutation properties 
with D 

fD + D7 5 = 0, 7 5 £> + Df = 0, (28) 



VD-DV = 0, OD-DO = 0. (29) 
This allows to define two sets of chiral projector operators 



P± = \(l±f), P ± = ^(l±7 5 )- (30) 

Similar pairs of projectors play a central role in approaches to chiral 
gauge theory based on the GW relation (@, @) and it's generalizations (|§) 



as domain wall fermions [0, Neuberger's overlap operator ||, proposals by 
Fujikawa @ and Chiu 0. 

On the other hand, it should be noted that our approach does not fit 
exactly into these schemes. The operators V, O do not have the spectral 
properties of the similar operators considered, for example, by |3J],|§. It 
seems that these differences may be understood as caused by different aims. 
The aim of the standard GW approach is to obtain a single Weyl fermion on 
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the lattice, without any doublers. In our approach we allow some doubling 
which gives a pair of Dirac fermions on the lattice, but we want nontrivial 
chiral symmetry to allow nontrivial chiral interactions between the doublers. 

If the properties of 7 5 ,7 5 we have found are sufficient to develop a con- 
sistent chiral gauge theory on this lattice remains to be shown. Note that 
the situation is further complicated by the non-pointwise character of the 
complex structure which remains to be understood. 
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